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We explore and calculate the rich scaling behavior of copolymer networks 
in solution by renormalization group methods. We establish a field theoretic 
description in terms of composite operators. Our 3rd order resummation of 
the spectrum of scaling dimensions brings about remarkable features: The 
special convexity properties of the spectra allow for a multifractal interpre- 
tation while preserving stability of the theory. This behavior could not be 
found for power of field operators of usual (j)^ field theory. The 2D limit of the 
mutually avoiding walk star apparently corresponds to results of a conformal 
Kac series. Such a classification seems not possible for the 2D limit of other 
copolymer stars. We furthermore provide a consistency check of two com- 
plementary renormalization schemes: epsilon expansion and renormalization 
at fixed dimension, calculating a large collection of independent exponents in 
both approaches. 

I. INTRODUCTION 

Recently, there has been considerable interest in the relation of field theory and mul- 
tifractals and the associated multifractal dimension spectra |]5|-|T0| as well as non- 



intersecting random walks and their 2D conformal theory |TT|-|I3|. We present a model of 
multi-component polymer networks that shows a common core of these topics and allows for 
a detailed study of the interrelations. While a multifractal spectrum may be derived from 
the scaling exponents of two mutually avoiding stars of random walks a description 

in terms of power of field operators seems to be ruled out by stability considerations 
Here we show that already a simple product of two power-of-field operators complies with 
both somewhat contrary requirements. We start from the theory of polymer stars, star - like 
arrangements of polymer chains with self and mutual excluded interaction [T^JT^. We gen- 
eralize this concept to stars of chains of different polymer species which may differ in their 
self and mutual interactions. Our formalism describes homogeneous polymer stars, stars of 
mutually avoiding walks (MAW), and the situation of two mutually interacting stars. 

Polymers and polymer solutions are among the most intensively studied objects in con- 
densed matter physics [1^. The behavior of multi-component solutions containing polymers 



of different species is especially rich. Systems of chemically linked polymer chains of dif- 
ferent species like block copolymers are also of considerable experimental and technical 
interest. Linking polymer chains of different or even contrary properties, like hydrophile 
and hydrophobe chains, one receives systems with qualitatively new behavior. Here, we 
concentrate on systems of chains in a solvent which, in a given temperature range, differ in 
their respective steric interaction properties. We generalize the theory of multi-component 
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polymer solutions |T9|-^ to solutions of copolymer networks, i.e. chains of different species 
linked at their endpoints in the form of stars or networks of any topology (see Figs.|l| - 
For solutions of polymer networks of a single species the scaling properties have been 



extensively studied by renormalization group methods (for a review see ||2^). Star polymers 
as the most simple polymer networks may be produced by linking together the endpoints of 
polymer chains at some core molecule (Fig. |l|). 

In the same way networks of any given topology may be generated (Fig. Q). Randomly 
linked polymer networks are also obtained as a result of a vulcanization process by randomly 
linking nearby monomers of different chains to each other. 

The asymptotic properties of homogeneous systems of linear chain molecules in solution 
are universal in the limit of long chains. Let us give a short account of the standard textbook 



results p!8| , p5| -p^] . For each system one finds a so called temperature at which the two 
point attractive and repulsive interactions between the different monomers compensate each 
other. As a result the polymer chains may be described by random walks (up to higher 
order corrections): The mean square distance between the chains endpoints (i?^) scales 
with the number of monomers N like (i?^) ~ A^. Above the B temperature the effective 
interaction between the monomers is repulsive resulting in a swelling of the polymer coil 
which is universal for N oo: 

{R") ~ N^-", (1) 

where the correlation length exponent u = 0.588, 3/4 in space dimension d = 3, 2. The 
number of configurations Z oi a. polymer chain of A^ monomers grows for A^ — > oo like 

Z ~ e^'^N^-^ (2) 

with a non - universal fugacity and a universal exponent 7 = 1.160, 43/32 for ci = 3, 2. 



In the early 70-ies following the work of de Gennes |^ the scaling theory of polymers was 
elaborated in detail using the analogy between the asymptotic properties of long polymer 
chains and the long distance correlations of a magnetic system in the vicinity of the 2nd 



order phase transition (see [[18| , P6| ). This mapping allows one to receive the above defined 
exponents u and 7 as limits of the correlation length exponent u and the magnetic suscepti- 
bility critical exponent 7 of the 0{m) - symmetric m- vector model in the formal limit m — 



On the other hand, if polymers of different species are present in the same solution the 
scaling behavior of the observables may be much more rich. Let us consider a solution of 
two different species of polymers in some solvent, a so called ternary solution. Depending 
on the temperature the system may then behave as if one or more of the inter- and intra- 
chain interactions vanish in the sense described above [p!9| -p3|. This will lead to asymptotic 
scaling laws that may differ from those observed for each species alone . 

Interesting new systems are obtained when linking together polymers of different species. 
The most simple system of this kind is a so called block copolymer consisting of two parts 
of different species. They are of some technical importance e.g. serving as surfactants 



3T| . For our study they give the most simple example of a polymer star consisting of 
chains of two different species (Fig. |^) which we will call here a copolymer star. For the 
homogeneous polymer star the asymptotic properties are uniquely determined by the number 
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of its constituting chains and the dimension of space . For the number of configurations 
(partition function) of a polymer star of / chains each consisting of monomers one 
finds fl^g: 

~ ei'^^Wf-^ ~ [Rjiy^s-i-^-' ,N^oo. (3) 

The second part shows scahng with the size R ~ A^'" of the isolated coil of monomers on 
some scale i, omitting the fugacity factor. The exponents ■yf,Tif, f = 1,2,3, .. . constitute 
families of 'star exponents', which depend on the number of arms / in a nontrivial way. 
The case of linear polymer chains is included in this family with the exponent 7 = 71 = 72 
defined in (0). For general numbers of arms / the star exponents 7/,?]/ have no physical 
counterparts in the set of exponents describing magnetic phase transitions. Nevertheless 
they can be related to the scaling dimensions of composite operators of traceless symmetry 
in the polymer limit m of the 0{m) symmetric m- vector model []32| , |33| . The exponents 
7/ have been calculated analytically in perturbation theory p6| , p!7|j3^ -p4| , by exact methods 



in two dimensions |1T7||35| , and by Monte Carlo simulations |36| , |37 |. 



It has been shown that the scaling properties of polymer networks of arbitrary but fixed 
topology are uniquely defined by its constituting stars ||2^ , as long as the statistical ensemble 



respects some conditions on the chain length distribution [^. Thus the knowledge of the 
set of star exponents 7/ or rjj allows to obtain the power laws corresponding to (|^) also for 
any polymer network of arbitrary topology. The partition function Zg of a polymer network 
Q (see Fig.|^) of F chains each of monomers scales with N —>■ 00 according to fTf^ : 

Zg ~ e'^^^AT^^^'e-^''^) , with rig = -dL + J2 ^fVf, (4) 

where n/ the number of vertices with / legs and L = l + I](//2 — is the number of 
loops in the network Q. 

In this article we address a somewhat more complex problem: What happens to the 
scaling laws if we build a polymer star or general network of chains of different species? In 
view of the above introduced ternary solutions, we thus study systems of polymer networks 
in which some of the intra and inter chain interactions may vanish. For instance, there may be 
only mutual excluded volume interactions between chains of the two different species in the 
copolymer star shown in |^ while chains of the same species may freely intersect. Each subset 
of chains of one species thus constitutes a star of random walks avoiding the second star of 
random walks. Gates and Witten [^] have shown that this problem of a star of random walks 
avoiding a given fractal structure is the key to calculating the the multifractal spectrum of 
absorption of diffusing particles on this fractal. This calculation can be performed explicitly 
for the diffusion near absorbing polymers |^,38 



The setup of our article is as follows. In section |T| we introduce notation and relate the 
polymer model to a Lagrangian field theory. This field theoretical formalism will be used 
throughout the paper. In section |IT| we define the renormalization group procedures. We 
present two alternative approaches: zero mass renormalization together with e- expansion 



(see e.g. [^]) and massive renormalization at fixed dimension |40|. Section IV is devoted 
to the study of the renormalization group flow of the ternary model and its fixed points. 
Series for critical exponents governing the scaling behavior of copolymer stars and stars of 
mutually avoiding walks are obtained in section 0. In section we discuss the problem of 
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resummation of the asymptotic series arising in this context. Numerical results are presented 
in section |VII| . We close with concluding remarks and an outlook on possible applications 
of the theory in section |V111| and give some calculational details in appendices. Some of our 
main results we have announced in a letter ET|. 



II. MODEL AND NOTATIONS 

Let us first take a look at the model we use to describe polymers. In a first discrete version 
we describe a configuration of the polymer by a set of positions of segment endpoints: 



Configuration{ri, . . . , r^v} G IR'^''^ ■ 

Its statistical weight (Boltzmann factor) with the Hamiltonian Ti divided by the product of 
Boltzmann constant k-Q and temperature T will be given by 

-1 1 ^ ^ 

exp[-— -K] = exp{--^ 5:(r, - r,_,f - ^ 5\r, - r,)}. (5) 

The first term describes the chain connectivity, the parameter io governs the mean sqare 
segment length. The second term describes the excluded volume interaction forbidding two 
segment end points to take the same position in space. The parameter (3 gives the strength 
of this interaction. The third parameter in our model is the chain length or number of 
segments N. The partition function Z is calculated as an integral over all configurations of 
the polymer divided by the system volume Q, thus dividing out identical configurations just 
translated in space: 

1 , .1 



= -J Y[dr,^xp[-^n{r,}]. (6) 

This will give us the 'number of configurations' of the polymer (|^). We will do our inves- 
tigations by mapping the polymer model to a renormalizable field theory making use of 
well developed formalisms (see [0,^ for example). To this end we introduce a continuous 
version of our model as proposed by Edwards |^,53] generalizing it to describe a set of 



/ polymer chains of varying composition possibly tied together at their end points. The 
configuration of one polymer is now given by a path r"(s) in d - dimensional space M!^ 
parametrized by a surface variable < s < Sa- We now allow for a symmetric matrix of 
excluded volume interactions Uab between chains a,b = 1, . . . , /. The Hamiltonian H is then 
given by 



,^^ni:n = i: t + 5 E I d^rp.(r)p.(r). (7) 



with densities pa{r) = Jq^" ds6'^{r — r"(s)) . In this formalism the partition function is 
calculated ClS db functional integral: 

Zf{S^} = I V[r%s)] exp{-^7^(r-)}. (8) 
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here the symbol P[ra(s)] includes normalization such that Z{Sa} = 1 for all Uab = 0. To 
make the exponential of 5-functions in @ and the functional integral well-defined in the 
bare theory a cutoff Sq is introduced such that all simultaneous integrals of any variables s 
and s' on the same chain are cut off by | s — s' |> Sq. Let us note here that the continuous 
chain model (|^) may be understood as a limit of discrete self-avoiding walks, when the length 
of each step is decreasing £o ^ while the number of steps Na is increasing keeping the 
'Gaussian surface' Sa = Nail fixed. The continuous chain model can be mapped onto 
a corresponding field theory by a Laplace transform in the Gaussian surface varables Sa to 
conjugate chemical potentials ("mass variables") /Iq PD| : 



Z}W}= / X{dS,e'^^'^Zf{Sa}. (9) 

The Laplace-transformed partition function Zf{^a} can be expressed as the m = limit 
of the functional integral over vector fields 0a, a = 1, . . . , / with m components 0^, a = 
1, . . . ,m : 

^fif^b} = J 'D[(f)a{r)] exp[-£{0fe, Hb}] \m=o- (10) 

The Landau-Ginzburg- Wilson Lagrangian C of f interacting fields (pb each with m compo- 
nents reads 



C{<Pi>,fii,} = \llj ^"^M + (V0.(r))2) + \Y. 

«=1 ■ a,a=l 



d''r<Pl{r)<Pl,{r), (11) 



here 0^ = Z]IILi(0a)^- The limit m = in (|TD|) can be understood as a selection rule for the 
diagrams contributing to the perturbation theory expansions which can be easily checked 
diagrammatically. A formal proof of ( |Tll) using the Stratonovich-Hubbard transformation 
to linearize terms in (|^) is given for the multi-component case in [^. The one particle 
irreducible vertex functions T^^\qi) of this theory are defined by: 

'^(E^OritLj?.) = / e^^'^Mri...dr^(0a,(ri)...0a,(r^))fpi,^=o- (12) 

The average (■ ■ ■) in ([T^) is understood with respect to the Lagrangian (^) keeping only 
those contributions which correspond to one-particle irreducible graphs and which have non- 
vanishing tensor factors in the limit m = 0. The partition function Z^f{Sa} of a polymer 
star consisting of / polymers of different species 1, . . . , / constrained to have a common end 
point is obtained from (|^) by introducing an appropriate product of (5-functions ensuring 
the "star-like" structure. It reads: 



Z.f{Sa} = J V[ra] exp{-^n{ra)} n S\ra{0) - n(0)). (13) 



n 

"B^ a=2 

The vertex part of its Laplace transformation may be defined by: 

S{P + E ?i)r^*^nP> ?i •••?/) = / e^(f"o+«»"'M'^rod'^ri . . . d'^rj 

(0i(ro) . . . 0/(ro)0i(ri) . . . 0/(r/))fpi^^=o, (14) 
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where all ai,...a/ are distinct. The vertex function T^*-f'> is thus defined by insertion of 
the composite operator Yia 'Pa- Its scaling properties are governed by the scaling dimension 
of this operator. When only one species is present one can also define T*^ by insertion 
of a composite operator of traceless symmetry In the following we will be mainly 

interested in the case of only two species of polymers, with interactions mh, U22 between 
the polymers of the same species and U12 = U21 between the polymers of different species. 
In this case the composite operator in ([T^) reduces to the product of two power-of-field 
operators with appropriate symmetry (0)-^^ (</>') ■^^ each corresponding to a product of fields 
of the same 'species'. Nevertheless, our results are easily generalized to the case of any 
number of polymer species. 

The starting point for our calculations are the three loop expansions for the bare vertex 
functions of interest [d / dk'^)V'^'^\ V^^\ Y'^*^\ They involve the loop integrals D2,Ii — Ig- 
These are given in appendix A together with their corresponding graphs. The expressions 
read: 

^ = 1 - l^2ul + ^hul, (15) 



8 , 40, 58, 14, 22, , 4 

^/3 + y/4 + ^/5 + ^/a + ^/vK„, (16) 

'^'^ 1 ~l~ ^2^ ^^0102/2 ~l" -D2ltaj^a2 ^(1304/8 ~l~ -1^122 ^0103-^1 ~l~ ^aiai^aia2-^l 
+ (Ji + D2) Maia2/2 + DlUa,a2Uaia4Ua5a6/^8 
~l"-D2ltaia2^aia3^a4a5-^l/2 ~l~ Ua^a2^a\az^a^a4^Ii 
~l~ (-^5 ~l~ -^g) ^0102^0103^0204/2 ~l~ Ua-ia2^a-iaz^a-ia4^I i 
+ (3 /4 + Jy) Uaia2 ^0103^^0203/3 + ^2 ijl + ^2) ^0i02^a3a4/4 
+ (D2/1 +2/4 + 2/5 + 17) Ul^a2^a,a, + {h + ^ h + h) ul^a,Ua,a2 
+ (D2h +2 h + h + h + Dl) Ul^jl + /^2MoioiMoi02Mo304A/2 
+ {Ji + -^5) "^0101*^0102 "^0203 + (-^4 + -^5 + -^7) "^0101*^0102 "^0103 

+MoiOiMoi02Mo2a2-^5/2 + (£'2/1 + 4 /4 + /e) Ua^a^ul^ai- i^^) 

In (|l^) summation over = 1 . . . / is assumed. The equations (|15|-|17[) apply to any number 
of polymer species. For a star of /i chains of species 1 and /2 chains of species 2 we restrict 
a = 1, 2 in (|15D,(|16D and the matrix of interactions Uab is given by 

r Mil I <a,b< fi 
Uil^' = \ U22 fi<a,b<f 
I ^12 else. 



Let us define in this way: 



Yi*hh) = Y(*f)\ .f,. (18) 

"ab— "jji, 
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For general /i, /2 the corresponding combinatorics may also be directly calculated by sum- 
mation over tti = 1,2 instead. Replacing Uaiaj = Ua^aj, each term in the sum with indices 
Oi . . . Ofc then acquires a factor 

f )f )• 

^#i(ai . . . Ofc)^ ^#2(01 • • • flfc)^ 

Here #i(ai . . . a^) is the number of Oj = 1 whereas #2(a'i • • • ctfc) is the number of = 2. 

As a special case we may derive the vertex function r^^22 the M12 interaction using 
the relation r*^*^^) = d / dui2TYi22 which is obvious from the perturbation theory (see [3C] for 
instance): 

rff22 = /dni2r(*^^). 

Note, that the vertex function T^*"^^^ and T^*^^^ define a vertex function with a 0^ insertion 
which in standard literature is denoted by F*^^'^^ |39|. With the same formalism we can also 



describe a star of / mutually avoiding walks [[T^jH • ^^^^ '^^^^ interactions on the same 
chain Uaa vanish and only those Uab with a ^ b remain: 

-pi*f) _ -pi*f)\ t-io\ 

■'-MAW — -'- \uab = {^-iab)ui2- \^^) 

In this case each term with indices ai . . . acquires a factor (^{^k\. 

As is well known, ultraviolet divergences occur when the vertex functions ( |15|) - (|17D are 



evaluated naively [Q. In the next section we apply the field theoretical renormalization 
group approach to remove the divergences and to make transparent the scaling symmetry 
of the problem. 



III. RENORMALIZATION 

We apply renormalization group (RG) theory to make use of the scaling symmetry of the 
system in the asymptotic limit to extract the universal content and at the same time remove 
divergences which occur for the evaluation of the bare functions in this limit p9| , P| , P5| . The 
theory given in terms of the initial bare variables is mapped to a renormalized theory. This 
is achieved by a controlled rearrangement of the series for the vertex functions. Several 
asymptotically equivalent procedures serve to this purpose. Here we will use two somewhat 
complementary approaches: zero mass renormalization (see for instance) with successive 
^-expansion |^ and the fixed dimension massive RG approach |^. The first approach 



is performed directly for the critical point . Results for critical exponents at physically 
interesting dimensions d = 2 and d = 3 are calculated in an e = 4 — c? expansion ||46|-49 



The second approach renormalizes off the critical limit but calculates the critical exponents 
directly in space dimensions d = 2, d = 3 |[50|J5T|] . It also gives quantitative results for the 



preasymptotic critical behavior |[5^ , |53[| . Most authors tend to prefer one method and to 
exclude the other for non obvious reasons. The application of both approaches will enable 
us in particular to check the consistency of approximations and the accuracy of the results 
obtained. 

Let us formulate the relations for a renormalized theory in terms of the corresponding 
renormalization conditions. Though they are different in principle for the two procedures. 
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we may formulate them simultaneously using the same expressions. Note that the polymer 
limit of zero component fields leads to essential simplification. Each field (pa, niass and 
coupling Uaa reuormalizes as if the other fields were absent. First we introduce renormalized 
couplings gab by: 

Uaa = fx''Z;j;^Zaagaa, a = 1, 2 (20) 
m2 = fi'Z^^'Z^^'Zugi2. (21) 

Here, /i is a scale parameter equal to the renormalized mass at which the massive scheme 
is evaluated and sets the scale of the external momenta in the massless scheme. The renor- 
malization factors Z^^, Zab are defined as power series in the renormalized couplings which 
fulfill the following RG conditions: 

d 

Z^S9aa)^'^fl{Uaa{gaa)) = 1, (22) 
Zfk^ {gaa ) ^aaaa {'^aa {gaa ) ) = /^^fl-aa, (23) 

Z^^{gn)Z^2{g22)T^n22{uab{gab)) = /i^5'i2- (24) 

These formulas are applied perturbatively while the corresponding loop integrals are evalu- 
ated for zero external momenta in the massive approach and for external momenta at the 
scale of /i in the massless approach as explained in appendix 0. In the massive case the RG 
condition for the vertex function F^^^ reads 

Z<pSgaa)'i^aa{Uaa{9aa))\k^=0 = a = 1, 2. (25) 



In the case of massless renormalization the corresponding condition reads [39 



Z^S9aa)'i^adiUaai9aa))\k^=0 = 0, a = 1, 2. (26) 

In order to renormalize the star vertex functions we introduce renormalization factors Z^j^ 

by 

4/"zi;:J'Z^f,j,T(*f^f^\Uab{gab)) = A-^^. (27) 

In the same way we define the appropriate renormalization for the vertex function of mutually 
avoiding walks (MAW): 

^(MAW/)rM^Aw(Wl2(^afe)) = • (28) 

The powers of n absorb the engineering dimensions of the bare vertex functions. These are 
given by 

Sf = f{e/2-l)+4- e. (29) 



The renormalized couplings gab defined by the relations (|20| ) ,(|21| ) depend on the scale pa- 
rameter /i. By their dependence on gab also the renormalization Z - factors implicitly depend 
on /i. This dependence is expressed by the RG functions defined by the following relations: 
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1^-7- Qah — Pab{ga'b')- 
Cl/i 


(30) 


Cl/i 


(31) 


/^^ln^*/i/2 = V*hf2i9ab)- 


(32) 


IJ,^ In Zmawj = vf^'^igab)- 


(33) 



The function t]^^ describes the pair correlation critical exponent, while the functions rj^^f^f^ 
and tff^ {gab) define the set of exponents for copolymer stars and stars of mutually avoiding 
walks. Note that Z^.2o renormalizes the vertex function with a (jP' insertion which coincides 
with r^*^°''. Consequently, the usually defined correlation length critical exponent u is ex- 
pressed in terms of functions 77^20 and rj^ (see next chapter). Explicit expressions for the (5 
and 7] functions will be given in the next section together with a study of the RG flow and 
the fixed points of the theory. 



IV. RENORMALIZATION GROUP FLOW AND THE FIXED POINTS: 
e-EXPANSION AND PSEUDO-e-EXPANSION 

Here, we want to discuss the RG fiow of the theory presented in section |T|. In particular, 
we want to find appropriate representations for the fixed points of the fiow. In a study 
devoted to ternary polymer solutions, the RG fiow has been calculated [^] within massless 
renormalization and is known to the third loop order of the e-expansion. Note, that for the 
diagonal coupling Qaa the corresponding expressions are also found in the polymer limit m = 
of the 0(m)-symmetric 0^ model. They are known in even higher orders of perturbation 
theory To third loop order the expressions read: 

= -^9aa + ^(4 + 25 + 2e^)gl^ -\(^ + t^e- llJe)gl, + 

1 



1/. 9x 9 1/21 215 

= -egaa + -(4 + 2e + 2e^)gl, - g (y 

^(79 - 22J + 33C(3))^7l + 0{gl), a = 1, 2. (34) 

I e 1 
PI,, = -egi2 + 3(1 + 2 + + ^22)^712 + -{2 + 6 + e^)gl^ - 

\{\ + ~ \j^){gli + ^22)^12 - ^(3 + - 3^£)(^ii + ^22)^?2 - 

1 1 1 27 

_(2 + 5£ - 2Je)gl^ + -(15 - J){gl, + gl^)g^2 + -{- + 

9C(3) - ^J){g\^ + gl^gl^ + - 3J)^7iito^??2 + ^(12 + 6C(3) - 

2J){gii + g22)gl2 + ^(6 + 3C(3) - 2J)gl^ + 0{g'>). (35) 

Here, the Riemann (^-function with C(3) ~ 1.202 and the constant J ~ 0.7494 occur. We use 
an index e at (3'^ to distinguish the /5-functions obtained in massless renormalization with 
successive e-expansion from (3'"^ obtained in massive field theory. 
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Similarly, performing renormalization in the massive scheme, we obtain the corresponding 
functions (3"^ . We present them using convenient variables Vab = D^Qab and introduce new 
functions = Here, is the one- loop integral calculated within massive field 

theory (see appendix A). This procedure defines a convenient numerical scale for the massive 
/^-functions. The expressions for the functions read: 

4?; 2 1 2 

= -(4 - d)vUl - ^ + g(22(zi - -) + 2^,)vl + -(-89 + 

310ii + 8^2 + Si2d - Uis - 180^4 - 87^5 - 21^6 - 33^7 - 

ms)vl) + 0{v'), a = 1,2. (36) 

-(^^11 +^^22 + 2t;i2) + -, 



(^m2 = -(4 - d)Vu{l - -{Vu + V22 + 2Vi2) + -(-^'ll - ^22 " ^v^Vn 



4f?2 o , 2^^?2«2 , 8vf^ii , , . , o 2 ■ , 2t;i2^ 



2f 12^^22 + —jf^ + ^^^^ + 4f i2t;22«i + 2vtiii + — ^ + 



3 - - 3 3 

Avuviiti + 2vl^%^) + ^ V'^'vi^vl^v'^^) + 0(t;^). (37) 

jfcZ 

Here, are the dimension-dependent loop integrals, normalized by the one-loop integral 
value (see appendix Expressions for the coefficients are given in the appendix 0. 
Note that the /3-function for the diagonal coupling gaa is known within the massive scheme 



to the order of six loops [E5| 



Let us solve the equations for the fixed points (FP) P({5'ii, (7225 5*12 }) of ^^e /5-functions, 

/?L(^?-) = 0, a = 1,2, 
/3,^,(^?n,^72*2,^7i*2) = 0. (38) 

As is well known, the first equation has two solutions (7*^ = 0,(73. -^^^ ^'^^ s^^" 
ond one finds a total of 8 FP depending on the choice of g*^^. The trivial FPs are 
Go(0, 0, 0), Uo(^g, 0, 0), U'o(0, gl, 0), So((?^, gl, 0), all corresponding to gl^ = 0. The non triv- 
ial FPs are found as G(0,0,^5), U(c/^, 0, ^{j), \J'{0,glg^), and S{glglg^). In the three 
dimensional space of couplings gu, 922, 912 these FPs are placed at the corners of a cube 
deformed in the gi2 direction (see figure Their e-expansions read | |30|| : 



sa = |-(^+5C(3))|^, (39) 
, 9c 39c^ I 267 693 189 n , 

,J.!£ + lM+(ilL-^C(3)-^j)e^ (41) 

4 128 V2048 256^ ^ 128 ^ ^ ^ 

For the evaluation of the fixed points of the /5-functions, calculated in the massive scheme 
(PB|), ( P7| ) (as well as of the other quantities of the theory) one has several alternatives. The 
first possibility is to introduce e-expansions for the loop integrals. For massive renormaliza- 
tion these are known for the one- and two- loop integrals (see ||45[] ): 



= ^(l - f ) + ^1 = ^ + 1 + 0{e\ ^2 = + 0{e\ (42) 
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A second possibility is to proceed directly in fixed dimension of space d = 2 or d = 3 
substituting numerical values for the loop integrals [^,^. For the massive two loop e- 
expansion the fixed points read: 

.s = x + ^, (43) 

* 3e 

.a = y + ^. (45) 

Only the first order of these results coincides with the fixed point values of the massless 
renormalization scheme |^ It is well known that the values of /5-functions, fixed points, 
and other intermediate functions in general depend on the RG scheme, only the critical 
exponents and other observables will be independent of the scheme followed. 

We will now study the expressions (p6|) , (^) directly at fixed dimension. In this scheme 
the usual way of finding the fixed points of /^-functions of models with several couplings 
involves the numerical solution of the system of equations (|3^) . To this end the asymptotic 
series in the coupling constants are represented in the form of corresponding resummed 
expressions (3'^'^^ |^ . However, the numerical solution of the resummed fixed point equation 
in general leads to inconsistent results, as we will show in Sec. [V^. An alternative to this 
procedure and thus a third possibility to proceed was originally proposed by Nickel and may 
be called a pseudo-e expansion |^ . As to our knowledge it has until now not been applied to 
theories with several couplings (see [^), although it seems a convenient tool to circumvent 
the specific difficulties arising for the massive approach. To apply this method, we introduce 
the "pseudo-epsilon" parameter r into the expressions for the /3-functions (p6D, 
(^) in the following way: 

-C7((4 - d)v^^) = r - + . . . , a = 1, 2, 

-/3,";^/((4 - d)v^2) = r- ^-{vn + V22 + 2^712) + • . . . (46) 

We solve for the fixed point solutions as series in r. The resulting series for the fixed points 
then either can be resummed (to obtain the numerical values of the fixed points) or they 
can be substituted into the expansions for the observables of the theory. In the final results 
we substitute r = 1. 

Performing this procedure we get the fixed point values as series in the pseudo-epsilon 
parameter r up to the order r^: 

9 9 ii , ^ , ^ ^ . 2 9 24 



VG = 3/2r + (3^1 - 3/2)r^ + (--^ _ ^ + 3/8 + 12 - - 

o 4 Z 

27^5 9^7 3 

r , 

8 8 ^ ' 

9r mil 93 3i2x 2 .387z6 9^8 33^2 128Hi 

Vu = h 1 T — 1 1 1 

8 ^ 32 64 32 ^ ^ 512 64 64 256 
459 693^7 27^3 369^4 1485 ^5 _ 3 ^2^ _ 9 Z2 ii _ 
512 ^ 512 ^ T28" ^ 64 ^ 512 64 8 



(47) 
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969 if 27di2 ,^3 



64 512 



)r^ (48) 

.Sio 33 zi 33, n , 27 261^5 

vs = 3 At + — + + 3/4 - 

* ^ U6 16 32^ ^ ^ 32 128 

135 u 261^1 63^6 3^2^ 33^2^1 9di2 99^7 9 is 

~32 64 l28" ^ ^ ^ 16 128" ~ l28" 

9^ 363^ 

32 32 ^ ^ ^ 

The expressions (^9]) - (^Tj) and the expressions (|36| ) - p9| ) give the fixed point values 
of ternary solutions in the massless and massive renormalization schemes and are the main 
results to be used in the subsequent calculations. 

Looking for the stability of the above described fixed points one finds that only the fixed 
point S is stable |^ . In the excluded volume limit of infinitely long chains the behavior of 
a system of two polymer species is thus described by the same scaling laws as a solution of 
only one polymer species. Nevertheless taking into account that real polymer chains are not 
infinitely long one may also find crossover phenomena which are governed by the unstable 
fixed points. Knowing the complete RG fiow allows to describe crossover phenomena in the 
whole accessible region However, for the purpose of our study we are interested only in 
the values of the fixed points and the properties of the star vertex functions at these fixed 
points. 



V. RESULTS FOR EXPONENTS 

For homogeneous stars of polymer chains of one species alone, several sets of star expo- 
nents have been defined, each describing either the scaling properties of the configurational 
number (see formula (^ of this article), or the anomalous dimensions of star vertices, etc. 



Due to scaling relations, these exponents can be expressed in terms of each other |2|]. In 
this sense, each set of star exponents forms a complete basis. For the copolymer and MAW 
stars, we here choose to present our results in terms of the exponents ^^/i/a and rjj^^^ given 
by the fixed point values of the functions r]^:f^f^{gab) (0) and rff^i^gab) (il)- Let us define 
the asymptotic values of copolymer star exponents and MAW star exponents by: 

V%f2=^*hf2{9ab)\s, (50) 

Vhf2='n*hf2i.gab)\G, (51) 

^/i/2 = ^*/i/2 {9ab)\u = V*f2h {9ab)\u', (52) 
^MAW^^MAW(^^^)|^_ (53) 

The exponent in the symmetrical FP S can also be expressed by f]fif2 — V^i+f2,o- Starting 
from the expressions for the fixed points given in the previous section, and the relations 
(^) - (|53|), we find the series for the star exponents. In e-expansion we obtain the following 
expansions for rjf-^f^: 

Vhd^) = -h + /i /2 (/2 - 3 + /i)^ - h /2(/2 - 3 + h)[h + 

/2 + 3C(3)-3)|^, (54) 
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^hh i^)=fi{^-fi-^f2)l+ fi (25 - 33/, + 8/, 2 - 91/2 + 42/,/, + 



18/2 



256 



+ fi 577 - 969/, + 456/,' - 64/, ^ - 2463/^ + 2290/,/, 



492/, 2/g + 1050/,' - 504/,// - 108// - 712C(3) + 936/,C(3) - 



224/, '((3) + 2652/gC(3) - 1188/,/gC(3) - 540//C(3) 



4096' 



(55) 



ifi-m^+fiifi 



l)(2/,-5)^-(/,-l)/,(4/,^ 



20/, +8/,C(3)-19C(3) + 25)-. 



(56) 



Here ({3) ~ 1.202 is the Riemann (^-function. The above formulas reproduce the 3rd order 
calculations of the scaling exponents of homogeneous polymer stars 7/ — 1 = i^iri^j^Q — fV2o) 
3^ . The exponents A*^^^-* given to 2nd order in equations (xx) of ^ to describe the multifrac- 



tal scaling properties of a Laplacian field with fractal boundary conditions are reproduced 
following A(29)(n) = -r^^n, A(^^)(n) = -r/^^^ + r/^o, Xi''Kn) = -r/^^, X(^^\n) = -r/U„, cor- 
recting a misprint in eq.(49) of 0. Also the 2nd order results for exponents x^^n — = 
-'^iVL,n - Vl,i) of [0] and the MAW exponents ctl = 1/2?]^^^ ^jefined in ^ find their 3rd 
order extension by the above expansions. The pseudo-e expansions for fjfif^ obtained in the 
massive scheme read: 



Vhh = ^ { - 1 + (/2 - 3 + /i) [r (zi - 1/2) + ^ (32 (zi - 1/2)2 + g ,g _ 18 _ 8 + 



22^i-627-6^5-(/2-3 + /i)(2 + 6«4-6«i))] - -(/1/2 - 2)(1 + 3 ^5 + 3 - 6 zi)}, 

(57) 



^{128 - 128 /i - 384 /a + r[(288 h - 144) /s' - 208 + 416 ii + 32 ^2 + 



(272 - 32 ^2 - 544 ii) f\ + (-32 ^2 - 1472 ii + 736) /2 + (-64 + 128 ii) // + 
(-336 + 672 zO hh] + ^ fi' f^'VuM-M)^ 

k\k2 



(58) 



^MAW 



T e 



/( /-I) 
16 



{ - 4 - (20 ii - 10 + 4 / - 8 fii)T + [(-3 + 18 ii - 12 ^4 - 3 



5)/^ + (3 Z5 + 32 i/ - 8 + 15 - 88 ii + 42 u + 23)/ - 26 - 18 zq + 19 k + 12 is + 



106^1 - 80 i/ - 30 ^4]^^}- 



(59) 



The expressions for the three loop terms rju-ki-M in (H) are given in the Appendix B. It has 
been pointed out in p that for the exponent 77^ = — A*^^^''(l) (see above) an exact estimate 
equal to our first order contribution may be found. It is indeed remarkable that all higher 
order contributions to appear to vanish in both approaches. 
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With these exponents we can describe the scahng behavior of polymer stars and networks 
of two components, generahzing the relation for single component networks In the 

notation of (^) we find for the number of configurations of a network Q of Fi and F2 chains 
of species 1 and 2 

Zgr^{R/ty^~^^^-^-^-^^\ withr^g = -rfL+ Nf,f,VM., (60) 

/l+/2>l 

where L is the number of Loops and Nf-^f^ the number of vertices with /i and /2 arms of 
species 1 and 2 in the network Q. To receive an appropriate scaling law we assume the 
network to be built of chains which for both species will have a coil radius R when isolated. 

For the sake of completeness we give also the results in e- and pseudo-e expansions for 
the correlation length critical exponent u = l/{2 + r]2o — f](j>) and the pair correlation function 
critical exponent rj^ in the non trivial fixed point: 

rije) = — + , (61) 

' 64 1024' ^ ^ 

/X 1 1 15 2 / 135 33C(3)\ , 

^ ^ 2 16 512 V8192 1024 y ' ^ ^ 

M't) = ^%^^(16r«s - 12x^^2 - 2AtHs + Sr^z^^ + SSr^z^zg), (63) 

, , 1 4- d (4:- d) , „ o 

z/ r = - + r - ^ ^ 4 - 40ii + 8^2 + 4d - 

^ ^ 2 16 512 ^ ' 

^ \ Q + 12^2 + lOii + 84z4 + 63^5 + 3z6 + 33«7 - 12^8 - 

220^1 + 4^2 + 2421^2 - + lOiid - 2i2d - y (64) 



VI. RESUMMATION 

As is well known, the perturbation series expansions of renormalized field theory are non- 
convergent but generally assumed to be asymptotic. For the exponents 77/^/3 this behavior is 
indicated in the corresponding columns of the tables | and 0, where the series are summed 
without further analysis. 

The increase in the coefficients of the high order terms of perturbation theory series may 
be estimated using information such as the combinatorial growth of the contributions with 
order. The series for the f3 function of the 0{m) symmetric 0^ model with one coupling g 
has the following asymptotic behavior 

(3{g) = J2A,g', (65) 

k 

Ak = ck'"'i-a)''k\[l+0{l/k)] ,k^oo. (66) 

The quantities a, Bq, c were calculated in pD|,B2|. A similar behavior is found for the critical 
exponents expressed as a series in powers of the coupling. The same results also apply to 
the divergence of the e and pseudo e expansions derived above. The property (B^) indicates 
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the Borel summability of the series f3{g) The Borel resummation procedure takes into 
account the asymptotic behavior of the coefficients and maps the asymptotic series to a 
convergent series with the same asymptotic hmit. The function I3aa (0) coincides with the 
0{m) symmetric (3 function (^) in the polymer hmit m = 0. So its asymptotic behavior is 
known. The asymptotic behavior of the off diagonal (3 function jSu was found by instanton 
analysis (see [|64| , |65| ) in []30| . 



Let us introduce the techniques for resummation, using the known asymptotic behavior 
of the series. Here, we apply Pade-Borel resummation |5y] and a resummation refined by 
a conformal mapping [^. The first way of resummation is applicable only for alternating 
series, while the second one is more universal. The resummation procedures are as follows 
50| , |5llj63| , |66[| . For an asymptotic series 



/(^) = E/^'V, (67) 
j 

one defines the Borel-Leroy transform f^{s) of the series by: 

with the Euler F-function (6 is a fit parameter). Then, the initial series may be regained 
from 

r%e) = / dtt'e-'fiet). (69) 

Substituting for f^{e) its analytic continuation in form of a Fade approximant and evalu- 
ating (^) for the truncated series, this procedure constitutes the Fade Borel resummation 
|50| , |63[] . The conformal mapping technique in addition uses the constant a in (|66D. Assuming 



the behavior (66) holds also for the expansion of f{e) in e, one concludes that the singu- 
larity of the transformed series f^{e) closest to the origin is located at the point (— 1/a). 
Conformally mapping the e plane onto a disk of radius 1 while leaving the origin invariant, 

_ (1 + aey/^ - 1 _ 4 w 

substituting this into f^{s) and expanding in w we receive a series defined on the disk with 
radius 1 on the w plane. This series is then resubstituted into (p9|). In order to weaken a 
possible singularity in the w-plane the corresponding expression is multiplied by (1 — w;)" 
introducing an additional parameter a In the resummation procedure the value of a 
is taken from the known high-order behavior of the e-expansion series while a is chosen in 
our calculations as a fit parameter defined by the condition of minimal difference between 
resummed 2nd order and 3rd order results. The resummation procedure was seen to be 



quite insensitive to the parameter h introduced in the Borel-Leroy transformation ( |68D [|5T 

For the resummation of the exponents Tjf^j^ we take into account the combinatorial factors 
which multiply each contribution according to the numbers of chains /i and f2- We include 
an additional factor (/i + f2)^ for the fcth order contributions, multiplying the constant a by 
(/i + /2). For resummation of the series at the fixed points S,G and U the following values 
oi a = ,a'^ ,aP are used [|60| , p0|] : 
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a^ = a^ = 3/8 and a" = 27/64. (70) 

By analogy we use the same procedures developed for the e expansion also for the r- 
expansion which we assume to have the same asymptotic behavior as it is in the same 
way collecting contributions of the same loop order. 

Let us note here that the conventional resummation of the (3- functions in the massive 
approach leads to a severe inconsistency which is the reason for us to take the pseudo-e or 
r-expansion method. The distinct feature of the /9a6-f unctions introduced here is that they 
are functions of different numbers of variables which leads to ambiguities in their analytical 
continuation via Pade approximants or rational approximants of several variables (see [|67[] ). 
Let us illustrate this for the example of the two-loop approximation. The corresponding 
expressions read: 

Pv^a = -(4 - d)VaaU^iVaa), a = 1, 2, (71) 
Pvi2 = -(4 - d)Vi2fvi2{Vu,V22,Vi2), (72) 

with obvious expressions for /i,^^ , f^^^ , /^^j • order to obtain the analytical continuation of 
the Borel transformed functions of one variable fvaai'^aa) (|7lD one can make use of the [1/1] 
Pade approximant. Solving the corresponding non linear equations numerically we find for 



the non-trivial fixed point S: Vu = V22 = 1.1857 ^8|. In order to apply a similar resumma- 
tion technique to the function /i,i2('yii, "^22, "^12) ( [T2| ) one can make use of a generalization of 
Pade approximants to the case of several variables, i.e. represent the Borel transform 



of /^j2 in the form of a rational approximant / of three variables 

rp f . . ,N 1 + aiivn,V22,Vi2)t + a2ivn,V22,Vu)t'^ 

[Vut, V22t, Vut) = — TT • (73) 

^ + b{Vn,V22,Vi2)t 



In spite of the fact that the rational approximant ( [73D preserves the projection properties of 
the initial series (0), i.e. putting any pair of variables {fn, ^22, ^^12} to zero in (|73D one gets 
the appropriate [1/1] Pade approximant for the remaining variable, the "global" symmetry 
is not preserved. Due to different analytical continuations for the Borel transforms of the 
series (^) on one hand and of (^) on the other, solving the fixed-point equation for the 
resummed function 

/;t = (74) 

we will never obtain a symmetrical solution vli = V22 = ^^12 7^ 0- For the fixed point 5" 
we substitute vn = V22 = 1.1857 and solving ([T^) we receive V12 = 0.9571 (!) [^. The 
reason is, that substituting numenca/ values of fixed point coordinates vu,V22 into (^) we 
lose information about the contributions to the fixed point value from different orders of 
the perturbation theory series. So it appears quite natural to restore this information by 



generalizing the pseudo-e expansion ||59| to the case of several couplings as described in 
section 0. 



VII. NUMERICAL RESULTS 



In the following we present our numerical results for the exponents ^^/j, ^X/a rj¥^^- 
The exponent in the symmetrical fixed point 5* is included due to the relation 77^^^^ = Vf^+fifl 
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71| . Numerical results for the exponent 7/ — 1 = I'ifjfQ — f'r]2o) be found in e-expansion 



m 



32| and in pseudo-e expansion in ||72 



A. d = 3 

Let us first consider the case d = 3. The tables || and |I| show some of the resummed 
results for the e- and r expansions in comparison with the naive resummation of the series. 
While the non-resummed results differ to a great extent for the two approaches at high /i, /2, 
resummation shows that the two schemes yield consistent numerical estimates. Tables |T|, 
IV| , and list our final results using the resummation procedure refined by the conformal 



mapping technique as described in the previous section. 

Comparing the numerical values listed in the above tables it is convincing that the two 
approaches and the different resummation procedures all lead to results which lie within a 
bandwidth of consistency, which is broadening for larger values of number of chains f\, f2 > 
1. This is not surprising as we have seen in section that our expansion parameters are 
multiplied by /i and f2- Rather it is remarkable that even for a total number of chains of the 
order of 10 (see tables p| , |rv| ) we still receive results which are comparable to each other. 



It seems noteworthy that at least for low numbers of chains (/i + /2 ~ 4) the non re- 
summed r-expansion seems to give results which do not differ essentially from the resummed 
values. Also the non refined Pade-Borel results of the r-expansion are closer to the refined 
summation of the e-expansion. 

Does the data answer the question of convexity of the spectrum? A close study of the 
matrix of values reveals, that for fixed fi both r^^^^ ^'^d r/y^j^ convex from above as 
function of /2, thus yielding 'MF statistics'. The relation to a MF spectral function for 
/i = 1, 2 has been pointed out in 0, it is analyzed in close detail in view of the new data 
and FT formulation in a separate publication ||3^. On the other hand also copolymer stars 



should repel each other. This is found to be true as well, the corresponding convexity from 
below shows up e.g. along the diagonal values rjff as function of /. The general relation 
^/i/2 + Vf[f2 — Vfi+f[j2+f2 always fulfilled. Thus, even though simple power k of field 
operators 0^ do not describe MF moments |jl[], they may be written as a power L + k of field 
operators of suitable symmetry which have the appropriate short distance behavior. This is 
also illustrated in the next subsection by Fig. ||, showing the spectrum of exponents ilfif2 
the 2D limit. The opposite convexity along the diagonal as opposed to each of the two axes 
is clearly seen for these combinations of two random walk stars which mutually interact. 



B. d = 2 

While two dimensional star polymers up to now have not found an experimental re- 
alization, their study is of some theoretical interest. It has been shown that the scaling 
dimensions of two-dimensional uniform polymer stars belong to a limiting case of the so- 
called conformal Kac table ||73|-[75|. They have been calculated exactly by Coulomb gas 



techniques [0,^. An exact relation has also been proposed for stars of mutually avoiding 
walks |Tl],IT2|. But it is still an open question if exact results for the copolymer star system 
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may be derived in this formahsm. Our numerical results for the exponents ^7/^/2, ''^/i/j; 
ffj^f^ are presented in tables 0, [VI I| , and 0. 



Exact results for exponents of two dimensional systems which are described by a con- 
formal field theory with central charge c < 1 may be taken from the Kac table of scaling 
dimensions ||7^-^: 

[(m + l)p - mqf - 1 
4m[m + 1) 

where p, q are integers in the minimal block 

1 < p < m — 1, I < q < p, (76) 
and m is connected with the central charge c by 

c = 1 — 6/m(m + 1), m > 3. (77) 
The exact result for the star exponents of uniform stars in two dimensions is received in the 



sub- limiting case of m = 2 (which means c = 0) for half integer values of p [|l7| , p5| : 

xj = 2/1^/2,0 = (9/2 -4)/48. (78) 
The scaling dimension x/ is related to the exponent rjf by: 

Xf = ^f{d -2 + ^)-^f. (79) 
For the exponents of the star of MAW the following result was conjectured for d = 2 



11,12.14 : 



^MAW ^ ^MAW ^ 2/,o,, = (80) 



These values are shown in the last column of table 0. Plotting the resummed data for rf^^ 
from table with respect to p one finds good agreement with the conjectured slope of 
-1/3. 

The qualitative behavior of the exponent Vjf^^^ in the Gaussian fixed point is shown in 
Fig. ^. The steps in the 'flying carpet' correspond to the difference of the results of the 
two RG approaches. Note that the curvature of the surface along the diagonal in the /i, /2 
plane has opposite sign to that along each of the axes. From this curvature it is obvious that 
the dependence of the exponent on /i,/2 may not be described by a simple second order 
polynomial. The best fit we could find to our resummed data using a simple formula which 
reproduces the vanishing result for /i + /2 = 3 found in e-expansion reads: 

^?;r = -/i/2(« + V(/i + /2)), a = 1/4, 6 = 3(1 -a). (81) 

Note that the right hand side of (|81|) vanishes if /i or /2 is zero according to our perturbative 
results. This might be a defect of the perturbation theory as a finite result may be expected 
in (i = 2 as in equations (|78D , (|80|) evaluated for / = 0. 

In 2D however, each chain of a star will interact only with its direct neighbors. A 
star described here by r/j^ will behave like a MAW 2/-star if each species-1 chain has two 
neighbors of species-2 whereas it will behave differently if the chains are ordered such that 
each species is in one bulk of chains. The 2D copolymer stars in this sense reveal an even 
richer behavior. Thus, the copolymer generalization of the MAW star adds another problem, 
for which a rigorous formulation in terms of an exactly solvable 2D model is yet to be found. 
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VIII. CONCLUSIONS AND OUTLOOK 



Several reasons motivated our study. First, we intended to reveal the scaling behavior 
of copolymer stars and networks in solutions generalizing former studies of homogeneous 
polymer networks. This included revisiting the theory of ternary polymer solutions and 
adding an independent approach to the calculations. Secondly, the description of multi- 
fractal spectra in terms of random walks promised to prove the relation of field theory 
and multifractals for this case. In particular we intended to check the convexity properties 
expected for the spectrum of exponents. A third motivation arose from the known pecu- 
liarities of polymers and polymer stars in two dimensions. Apart from numerically verifying 
previous results on polymer and mutually avoiding walk stars we pose the question of finding 
an exactly solvable (conformal) two dimensional theory for general copolymer stars. 

We have extensively studied the spectrum of exponents governing the scaling properties 
of stars of walks taking into account the self and mutual interactions of a system of species 
of polymers. Our study was performed in the framework of field theoretical RG theory 
using two complementary approaches: The renormalization at zero mass in conjunction 
with e expansion and massive renormalization at fixed dimension with numerical evaluation 
of loop integrals. We have formulated the problem of finding the scaling exponents of stars of 
interacting and non interacting walks in terms of the determination of the scaling dimensions 
of composite field operators of Lagrangian field theory. On the one hand this allows for the 
application of well developed formalisms and methods for analyzing the scaling properties. 
On the other hand this defines these new families of exponents extending previous sets in the 
framework of Lagrangian field theory. Our results agree with the previous studies of special 
cases which were in part done only to 2nd order of the e expansion. We have here considered 
the general case of a star of two mutually avoiding sets of walks, the walks of each set either 
self- interacting or not. Also we have studied the case of a star of mutually interacting walks. 
All calculations were performed to third order of perturbation theory. The sets of exponents 
are given in e-expansion (formulas ( [5^ ) - (|56|)) and in terms of a pseudo-e-expansion ( |57[ ) 
- (59)). The latter has proven to be a most suitable tool to evaluate this massive theory 
containing serveral couplings. We have shown that the conventional way of direct solution 
even of the resummed expressions for the fixed points of the theory would lead to severe 
problems in this case. We have evaluated the series obtained in both approaches for space 
dimensionality d = 2 and d = 3. Numerical values are produced by careful resummation of 
the asymptotic series using the results of an instanton analysis of the three coupling problem 
||3CI|| . For comparison we have also given the results of naive summation as well as standard 
Fade Borel resummation for selected cases. 

We have found remarkable consistency and stability of the results in ci = 2 and d = 3 with 
expected growing of deviations for large number of arms of one star. The same methods 
were applied previously to the problem of uniform star polymers and have led to results 
||3^ , [72[| in good agreement with Monte Carlo (MC) simulations |]36|j37[] . We hope our present 
calculations might also stimulate MC studies of the copolymer star problem. 

The study we performed for two dimensions might have no direct application to the 
physics of real polymers but it could perhaps give some insight to the problem of mapping our 
theory to a two dimensional conformal field theory. The resummed values of the exponents 
for stars of mutually avoiding walks are in fair agreement with an exact result previously 
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conjectured [|TT],|T2|. The exponents for the case of stars of two mutually avoiding sets of 
walks on the other hand show a dependence on the numbers of walks which may not be 
described by a second order polynomial as derived from the general Kac formula |73|-|75 



This may be seen already qualitatively from the fact that the curvature of the function ^7/1/2 
(see Fig. ^ of the two variables /i, /2 along each of the axes in the /1/2- plane has the 
opposite sign to the curvature along the diagonal fi = f2- 

It is this fact on the other hand, which shows that the series of exponents ^7/1/2 is a good 
candidate for finding its application in the theory of multifractal (MF) spectra ||7^ . The MF 
spectrum describing the moments of a fractal probability measure fulfills exact conditions 
of convexity. Deriving such a MF spectrum on the other hand from the scaling dimensions 
of a series of composite field operators is only possible if the scaling dimensions show the 
appropriate convexity |]T|. This in fact is given for our case and the series of exponents may 
be related to the MF spectrum generated by harmonic diffusion near an absorbing fractal 
1^. This also allows for a field theoretic test of results for the short-distance correlations on 

in 



multifractals 



This relation and the calculation of the MF spectrum on the basis of the 



results presented here, is subject of a separate publication pH . 
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APPENDIX A: LOOP INTEGRALS. GRAPHS, e-EXPANSION, NUMERICAL 

VALUES 

This appendix is devoted to the contributions to perturbation theory, their representation 
in terms of Feynman graphs and their corresponding loop integrals, and the evaluation of 
these integrals for the two RG approaches. Fig.^ shows the Feynman graphs up to third 
loop order representing the contributions to the functions F*^^) and F*^^) (we keep labeling of 



Each contribution to F*^*-^-' contains the composite operator nf=i 4'ai only once. The 
relevant graphs can be obtained from the usual four-point graphs 2-U2 - 12- U4 by considering 
each vertex in turn to describe the composite operator. In the three-loop approximation we 
consider here two more graphs contribute which can not be produced in this manner. They 
are labelled as 13 and 14. In table |VI11| we show the correspondence between the numerical 



values of the loop integrals and appropriate Feynman graphs. 
A diagram with L loops is to be multiplied by with 



2d-i(7r)d/2r(rf/2)' 

but this factor can be absorbed by redefinition of the coupling constants gab gab/ era- 
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In the massless renormalization scheme loop integrals corresponding to these graphs 
are evaluated by e-expansion at zero mass and non-zero external momenta chosen at the 



so-called symmetry point. The expressions read [p9|] : 



Dl = 


l/^ e e'^\ _ 1 3e be 
:-! + - + -, /^^ = — 1 + — + — 
2 2^ 2£:2V 2 2 


2 ^ 


11 = 


1 be\ „ 1 15£:\ 




11 = 


1 / 3l£2 3Je% 

6.3(1+3.+ 4 - 2 ). 




11 = 


1 he 23^2 3J£2^ 1 I 
3£3V^+ 2+4 2 ^ 6"3£3i 


1.2..-/) 


!',= 


25' 6e2l^ + ^^J- 





(A.l) 

Here, the values of derivatives d/dk'^ of the function r(^^(/c) are given at the point k'^ = 1. 
In the massive renormalization scheme loop integrals are calculated at non-zero mass and 
zero external momenta (to distinguish from ( [A.l|) we will label them by "m"). The mass 



renormalization introduces a higher order correction to the propagator, which has to be 
taken into account in our calculation only in the first order term (see eqs. (|16D, (p!?!)): 



Here D21 = (4 — rf) /4. This value has been substituted into the results for the beta functions 
and fixed points. D21 does not enter expressions which are independent of the RG scheme, 
such as the resulting exponents. The integrals can be either e-expanded (see formulas ( ^2[ ) 
from this article for instance) or numerically calculated at arbitrary space dimensions |]56| , |57| . 
In particular, for dimension d = 2 and d = 3 they are given in table ^ with the following 
normalization: 

Note that in the massive scheme the values of the derivative d/dk^ of the function r(^)(/c) 
are given at fc^ = 0. 



APPENDIX B: THREE LOOP CONTRIBUTIONS 

In this appendix we have collected the more lengthy expressions for the three loop con- 
tributions to RG functions and exponents. The coefficients (j -|- -|- / = 3) for the r 



expansion of the function /S^^ (pTQ read: 

^300 ^ ^030 ^ 



12 d 2^3 1624 iq 4^8 58 ii 2x2 20 



18 9 9 3 9 27 3 27' 



003 _ _824 _ ^ _ ^ _ ^ , 34 Zi 

~ 9 9 9 27 27 ' 

6210 = 5120 ^ Q 
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,102 ^ ,012 ^ 20 u _ ii? _ 2|6 _ 26 70 zi 
~ 9 9 9 27 27 ' 

L201 _ 1^021 _ _ 16^4 _ ^ _ ^ _ ^ 20^1 2^2 

~ ~ 9 9 3 18 27 9 27 ' 
,111 _ _ 16 2-^6 8z4 Ail 
~ 27 9 9 3 ■ 

The coefficients 'r]u;ki;k2 introduced for the r expansion of the exponent v'fifi unsym- 
metric fixed point U (^) read: 

r]u-ofi = 328 - 1480 ii - 128 is - 240 u - 492 ig + 132 ig - 356 h - 48 ig + 2288 ii^ + 
16^2^ + 384^1^2, 

W;0;i = -7680 ii^ - 496 ^1^2 + 5708 ii + 184 Z2 + 1326 h + 1620 - 588 iq + 204 h - 
16^2^ + 48^8 - 1312, 

riu.fl,2 = 570 + 810 u + 1488 ii^ - 24 ^2 + 48 ^1^2 - 2154 h - 270 ir + 216 ie - 216 25, 
^?c/;0;3 = -54 - 162^4 + 162^1, 

r^^.-^.Q = -16 ^2^ + 468 ij - 448 ^1^2 + 160 ^2 + 2408 h - 252 - 2992 ii^ + 564 + 
48^8 - 560, 

r^^.-^.^ = -5342 ii + 176iii2 - 88*2 + 756^6 - 594^7 - 252^5 + 3536 ii^ + 1638*4 + 1346, 
r]u;i:2 = 1188 ii - 216 - 216 - 756 u - 324, 

Vu;2;o = -1072 ii - 32 ^2 + 64 ^1^2 + 336 ^4 - 112 ^7 + 144 + 704 h'^ - 48 + 272, 
Vu;2;i = 1098^1 - 180*6 " 738^4 - 180^5 - 306, 
ryj;.3.o = -40 + 144 ii - 24 is - 96 u - 24 zg. 
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TABLE 1. Values of the copolymer star exponent rjf^f^ obtained in the first, second, and third 
order in the Gaussian (G) fixed point in e-cxpansion and pseudo-e (r) expansion for different values 
of /i,/2 at e = 1 {d = 3). res stands for the results obtained by Pade-Borel resummation of the 
three loop series. 



fl 


/2 


^ £ 






res. 


~ r 






res. 


1 


1 


-.50 


-.63 


-.46 


-.56 


-.50 


-.58 


-.56 


-.57 


1 


2 


-1.00 


-1.00 


-1.00 




-1.00 


-1.00 


-1.00 




1 


3 


-L50 


-1.13 


-1.99 


-1.36 


-1.50 


-1.25 


-1.42 


-1.34 


2 


2 


-2.00 


-1.50 


-2.65 


-1.81 


-2.00 


-1.67 


-1.93 


-1.80 


2 


3 


-3.00 


-1.50 


-5.71 


-2.50 


-3.00 


-2.00 


-3.01 


-2.45 


3 


3 


-4.50 


-1.13 


-12.27 


-3.48 


-4.50 


-2.25 


-5.09 


-3.37 


TABLE IL 


Values of the copolymer star exponent ///^/a 


obtained in 


the first, 


second, ; 


and third 


order 


in the unsymmetrical (U) fixed point in e- 


-expansion and pseudo-£ 


: (r) expansion for different 


values of /i,/2 


at e = 


1 {d = 3). 


res stands for the results obtained by Pade-Borel resummation 


of the 


ilircc loop Hvrici 
















fl 


h 


^ e 






res. 


~ r 






res. 


1 


1 


-.38 


-.50 


-.28 


-.43 


-.38 


-.46 


-.43 


-.44 


1 


2 


-.75 


-.85 


-.69 


-.80 


-.75 


-.82 


-.78 


-.80 


1 


3 


-1.13 


-1.07 


-1.33 


-1.11 


-1.13 


-1.09 


-1.11 


-1.10 


2 


1 


-1.00 


-.98 


-.71 


-1.00 


-1.00 


-.99 


-.98 


-.99 


2 


2 


-1.75 


-1.37 


-2.37 


-1.62 


-1.75 


-1.50 


-1.71 


-1.60 


2 


3 


-2.50 


-1.47 


-4.99 


-2.19 


-2.50 


-1.82 


-2.56 


-2.13 


3 


1 


-1.88 


-1.28 


-1.70 


-1.50 


-1.88 


-1.48 


-1.82 


-1.64 


3 


2 


-3.00 


-1.36 


-6.19 


-2.47 


-3.00 


-1.91 


-3.18 


-2.43 


3 


3 


-4.13 


-1.02 


-12.83 


-3.26 


-4.13 


-2.06 


-4.97 


-3.14 



TABLE III. Values of the copolymer star exponent vff^f^ at d = 3 obtained by £-expansion (e) 
and by fixed dimension technique (3d). 



/l 


1 


2 




3 




4 




5 




6 




/2 


e 3d 


e 


3d 


e 


3d 


e 


3d 


e 


3d 


e 


3d 


1 


-0.56 -0.58 


-1.00 


-1.00 


-1.33 


-1.35 


-1.63 


-1.69 


-1.88 


-1.98 


-2.10 


-2.24 


2 




-1.77 


-1.81 


-2.45 


-2.53 


-3.01 


-3.17 


-3.51 


-3.75 


-3.95 


-4.28 


3 








-3.38 


-3.57 


-4.21 


-4.50 


-4.94 


-5.36 


-5.62 


-6.15 


4 












-5.27 


-5.71 


-6.24 


-6.84 


-7.12 


-7.90 


5 
















-7.42 


-8.24 


-8.50 


-9.54 


6 




















-9.78 


-11.07 
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TABLE IV. Values of the copolymer star exponent ifj-^f2 at = 3 obtained by e-expansion (e) 
and by fixed dimension technique (3ci). 



/l 


1 




2 




3 




4 




5 




6 




/2 


e 


2,d 


e 


M 


£ 


3d 


£ 


3d 


£ 


3d 


£ 


3d 











-0.28 


-0.28 


-0.75 


-0.76 


-1.36 


-1.38 


-2.07 


-2.14 


-2.88 


-3.01 


1 


-0.43 


-0.45 


-0.98 


-0.98 


-1.64 


-1.67 


-2.39 


-2.47 


-3.21 


-3.38 


-4.11 


-4.40 


2 


-0.79 


-0.81 


-1.58 


-1.60 


-2.44 


-2.52 


-3.33 


-3.50 


-4.28 


-4.57 


-5.29 


-5.73 


3 


-1.09 


-1.09 


-2.13 


-2.19 


-3.16 


-3.30 


-4.20 


-4.48 


-5.28 


-5.71 


-6.41 


-7.03 


4 


-1.35 


-1.37 


-2.61 


-2.71 


-3.82 


-4.04 


-5.02 


-5.40 


-6.24 


-6.81 


-7.48 


-8.28 


5 


-1.60 


-1.64 


-3.05 


-3.21 


-4.44 


-4.75 


-5.80 


-6.30 


-7.15 


-7.89 


-8.51 


-9.50 


6 


-1.81 


-1.89 


-3.46 


-3.68 


-5.01 


-5.42 


-6.53 


-7.15 


-8.02 


-8.92 


-9.50 


-10.69 



TABLE V. Values of rjj exponents of star of mutually avoiding walks at d = 3, d = 2 
obtained by e-expansion (e) and by fixed dimension technique (3d, 2d). The last column gives the 
exact conjecture at d = 2 jri| , p^ ]. 





d = 3 






d = 2 




/ 


£ 


3d 


£ 


2d 


exact 


1 














-.250 


2 


-.56 


-.56 


-1.20 


-1.19 


-1.250 


3 


-1.38 


-1.36 


-2.71 


-2.60 


-2.916(6) 


4 


-2.36 


-2.34 


-4.36 


-4.07 


-5.250 


5 


-3.43 


-3.43 


-6.04 


-5.61 


-8.250 


6 


-4.58 


-4.64 


-7.78 


-7.17 


-11.916(6) 



TABLE VI. Values of the copolymer star exponent Vfif2 at d = 2 obtained by e-expansion (e) 
and by fixed dimension technique (2d). 

/i 1 2 3 4 5 6 



/2 


£ 2d 


£ 


2d 


£ 


2d 


£ 


2d 


£ 


2d 


£ 


2d 


1 


-1.20 -1.22 


-1.98 


-1.98 


-2.56 


-2.58 


-2.99 


-3.04 


-3.36 


-3.43 


-3.68 


-3.78 


2 




-3.41 


-3.45 


-4.49 


-4.59 


-5.37 


-5.52 


-6.13 


-6.34 


-6.80 


-7.04 


3 








-6.05 


-6.23 


-7.36 


-7.63 


-8.49 


-8.84 


-9.50 


-9.91 


4 












-9.06 


-9.44 


-10.55 


-11.03 


-11.89 


-12.45 


5 
















-12.38 


-12.98 


-14.03 


-14.74 


6 




















-15.99 


-16.81 
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TABLE VIL Values of the cc 


ipolymer star exponent t/^ at d 


= 2 obtained by e-expansion (e) 


and by fixed dimension technique (2d). 










/i 1 


2 




3 


4 




5 


6 


/2 e 2d 


e 


2d 


e 2d 


e 2d 


e 


2d 


e 2d 





-0.59 


-0.62 


-1.51 -1.53 


-2.61 -2.63 


-3.84 


-3.89 


-5.18 -5.28 


1 -.91 -.96 


-1.94 


-1.96 


-3.11 -3.13 


-4.35 -4.41 


-5.71 


-5.83 


-7.17 -7.35 


2 -1.62 -1.63 


-3.05 


-3.09 


-4.49 -4.54 


-5.94 -6.06 


-7.46 


-7.64 


-9.04 -9.30 


3 -2.16 -2.16 


-4.00 


-4.04 


-5.70 -5.80 


-7.39 -7.57 


-9.09 


-9.35 


-10.82 -11.17 


4 -2.60 -2.64 


-4.80 


-4.88 


-6.79 -6.96 


-8.72 -8.97 


-10.61 


-10.95 


-12.52 -12.95 


5 -3.00 -3.03 


-5.52 


-5.63 


-7.81 -8.01 


-9.96 -10.27 


-12.06 


-12.47 


-14.14 -14.65 


6 -3.34 -3.39 


-6.17 


-6.33 


-8.73 -8.99 ■ 


-11.12 -11.49 


-13.43 


-13.91 


-15.69 -16.27 


TABLE VIII. 


Values of the loop integrals. The graphs 2-Ml 


- 5-S3 indicate the derivative of 


the two-point function d / dk'^T^'^\k) . 










Graph 


Integral 




Graph 


Integral 




Graph 


J-ilLCCjl Oil 




value 






value 






value 


2-U2 


D2 




8-U4 


h 




14 


hD2 


3-U3 


Dl 




9-U4 


h 




2-Ml 





4-U3 


h 




10-U 


h 




3-S2 


h 


5-U4 






11-U4 


h 




4-M3 





6-U4 


I1D2 




12-U4 


I7 




5-S3 




7-U4 


h 




13 


Dl 








TABLE IX. Numerical values of normalized loop integrals ij 


calculated in ■ 


the massive field 


theory framework p6|,^. 














k 
















d = 2 0.781302412896 


-0.114635746230 -0.044703881514 


0.569829439192 


d = 3 0.6666666667 




-0.0740740741 


-0.0376820725 




0.3835760966 
















is 


d = 2 0.659043562065 


0.650899895132 0.40068563 




-0.157398409771 


d = 3 0.5194312413 




0.5000000000 


0.1739006107 




-0.0946514319 
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i^J ^0 ■ 

1 3 . ■ - 
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FIG. 1. Star polymer of / arms linked together at point ro with extremities placed at points 
ri...r/. 



c 



V \ y ^ 



FIG. 2. A polymer network Q. It is characterized by the numbers n/ of /-leg vertices. Here 
ni = 3, n3 = 2, = 1, ns = 1. 




FIG. 3. a. Block copolymer consisting of two polymer chains of different species (shown by 
solid and thin lines) linked at their endpoints. b. Copolymer star consisting of arms of species 
a and arms of species b tied together at their endpoints. 
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FIG. 4. Fixed points (FPs) of ternary polymer solution. The trivial FPs Go, Uq, Uq, Sq 
correspond to vanishing mutual interaction. The non-trivial FPs G, U, U', S correspond to 
non- vanishing mutual interaction (gi2 ^ 0). 




FIG. 5. Exponent ilf^^f^ the 'Gaussian' fixed point at c? = 2 obtained in e-expansion and in 
fixed d scheme. Steps on the "flying carpet" correspond to the difi'erence of the results of the two 
renormalization group approaches. The line across traces the diagonal values •q'^j. 



29 



1--50 2-Ml 3-52 4_7Vf3 5-53 1 - Ul 2 - U2 3 - US 4 - U3 

XxxXX>0>S( ^xfi'ysKun M 

5-C/4 6-C/4 7-C/4 8 - C/4 ^ _ V 10 - C/4 11 - C/4 ^3 - t/4 

13 14 

FIG. 6. Graphs of functions F^^^, F^^^ in three-loop approximation. Graphs 13 and 14 represent 
additional contribution to the function F^*-^). In 13, 14 the /-vertex is indicated by a box. 
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